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Observations of cosmic rays with energies above ~ 4 x 10^° GeV have 
^ ' inspired several speculative suggestions concerning their origin. The crucial 

I question is whether or not the spectrum exhibits the expected 'GZK cutoff' 

^SJ ■ at this energy — concerning which there are presently contradictory results. 

CNJ ' If there is indeed a cutoff, then the sources are cosmologically distant and 

O^l , rather exotic in nature. If there is no cutoff then new physics is required. 

m 

O [ 1. Introduction 

^ • As we approach the centennial of the discovery of cosmic radiation by 

Victor Hess in 1912, the astrophysical origin of these high energy parti- 
cles remains still unknown. Even so the study of cosmic rays has been 
, extremely rewarding for particle physics. As Hess noted in his Nobel Lec- 

ture in 1936: "The investigation of . . . cosmic rays . . . has led to the discov- 
. - ery of the positron . . . It is likely that further research into "showers" and 

I "bursts" of the cosmic rays may possibly lead to the discovery of still more 

• elementary particles ■ ■ ■ of which the existence has been postulated by some 

theoretical physicists in recent years". Indeed several such discoveries were 
made in rapid succession — the muon in 1937, the pion in 1947 and strange 
particles (kaons, hyperons) also in 1947. However with the development of 
particle accelerators in the 1950's, the attention of particle physicists turned 
naturally to controlled laboratory experiments. No other new particles have 
subsequently been discovered in cosmic rays although there have been sev- 
eral false alarms (e.g. free quarks, monopoles, . . . ) and many claims of 
unexplained phenomena (e.g. 'Centauro' events). Today we have another 
such mystery — the observation of ultra high energy cosmic rays (UHECRs) 
which cannot propagate very far through the cosmic microwave background, 
are unlikely to be significantly affected by cosmic magnetic fields, and yet 
which do not point back to any plausible astrophysical sources in our vicin- 
ity [I, 2]. The intriguing question is whether an explanation is possible in 
terms of known physics or will require physics beyond the Standard Model. 
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2. Observational Status 

As we have just heard from Prof Teshima, the AGASA collaboration has 
reevaluated the energy determination method used for analysing their data 
on air showers collected over 12 years, with particular attention to the lateral 
distribution, attenuation with zenith angle, shower front structure, delayed 
particles observed far from the shower core, etc [3]. They confirm that the 
energies assigned to AGASA events have an event-reconstruction accuracy 
of ±25% at 10^^ GeV, while the systematic uncertainty is ±18% (indepen- 
dent of primary energy above 10^'^ GeV). As seen in Fig. 1, the AGASA data 
at the highest energies connects smoothly with the (better sampled) Akeno 
data at lower energies. There are 59 events with E > Eqzk =^ 4 x 10^° GeV 
(and 8 events beyond 10^^ GeV) which have a spectrum distinctly flatter 
than at lower energies, suggestive of a different origin, but with no indica- 
tion of the GZK cutoff expected if the sources are cosmologically distant. 




log(ENERGY in eV) loq,c,(E) (eV) 

Fig. 1. The left panel shows the recalibrated AGASA spectrum along with data 
from other experiments at lower energies [3]. The right panel shows the HiRes 
data (fitted to a two-component source model incorporating a GZK cutoff for the 
extragalactic component); the AGASA data are shown for comparison, and with 
the event energies reduced by 20% [5]. (All fluxes have been multiplied by E^). 

By contrast the HIRES air fluorescence experiment, with a similar expo- 
sure in the mono mode, has reported only 1 event above 10^^ GeV, consistent 
with a GZK cutoff [4]. But as is clear from Fig. 1, the absolute fluxes are 
lower than those measured by AGASA and the begining of the 'ankle' in 
the spectrum is distinctly lower, suggestive of an energy calibration mis- 
match bewtween the fluorescence and air shower detector data. (Since the 
spectrum falls as ~ below the ankle, the fractional error in the flux 
is twice that in the energy, keeping in mind the change in the differential 
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energy interval with the energy.) If the AG AS A energies are lowered by 
~ 20%, the spectral shapes can be matched below 10^^ GeV [5], however 5 
AGASA events still remain above this energy. Given the low event statis- 
tics, the significance of the discrepancy between the two experiments is not 
overwhelming, nevertheless it is clearly of paramount importance to resolve 
the issue and establish a consistent energy scale. The engineering array 
of the Pierre Auger Observatory [6] has already collected several 'hybrid' 
events observed using both type of detectors so progress is expected soon. 

The other important new result concerns the AGASA data on the an- 
gular distribution of events on the sky [7]. As shown in Fig. 2, although 
the large-scale distribution of the 59 observed showers with E > Eqzk is 
consistent with isotropy, there are a number of 'clusters' — defined as a 
grouping of 2 or more events within (approximately the experimental angu- 
lar resolution of) 2.5*^. The chance probability that the 5 observed doublets 
and 1 triplet result from an isotropic distribution is estimated by Monte 
Carlo to be less than 10~^ [7]. However Fig. 2 shows that when the events 
are partitioned by energy, small-scale clustering is seen only for events with 
ii^ < 6 X 10^*^ GeV; at higher energies there is no clustering [8]. 




Fig. 2. Arrival directions of UHECRs of energy E > ix 10^" GeV (circles) and E > 
10^^ GeV (squares) , with the 'doublets' and 'triplet' highlighted (the shaded area is 
unobservable by AGASA) [7] . The bottom panel shows the angular autocorrelation 
function for events above various energy thresholds [8]. 



4 



To estimate the significance of tlie clustering reliably, one should calcu- 
late the 'penalty factor' for making a 'posteriori cuts (in energy thresholds 
and angular separations) on the data set in order to maximize the clustering 
signal. Secondly, as was emphasized earlier [9], the data set used to make 
the initial claim for clustering ought not to be used in the actual analysis. 
Recently both these issues have been investigated in detail [10]. As shown 
in Fig. 3, the strongest autocorrelation signal {P^^^ — 8.4 x 10~^) is seen 
at separation angle = 2.5*^, with energy threshold = 4.9 x 10^'^ GeV. 
However Monte Carlo tests show that 3475 out of 10^ simulated AGASA 
data sets have < P^i^ so the chance probability for this is 0.35%. 

This is 10 times higher than the value claimd earlier [11] — these authors 
did not allow for their (arbitrary) choice of 9c to maximise the clustering 
signal and interpreted their result as implying that the "Correlation func- 
tion of ultra-high energy cosmic rays favors point sources". Secondly, when 
the AGASA data is divided into 2 roughly equal sets (30 events detected 
before October 1995 and 27 events detected afterwards), the chance prob- 
abilities jump to 4.4% {0c = 2.4°) and 27% (61c = 4.7°) respectively for the 
observed clustering [10]! The reason appears to be that the two sub-sets 
of data are strongly correlated — the 'triplet' being made of 2 events from 
the first set and 1 from the second. Moreover AGASA has announced more 
events bringing to 42 the number detected since October 1995. There are 2 
'doublets' separated by less than 2.5° in this data set but the chance prob- 
ability for having observed these is 19% [10]. Thus the clustering observed 
by AGASA is not statistically significant, so does not require the observed 
UHECRs to originate from point sources. 

Nevertheless, some authors have sought for and found statistially signif- 
icant correlations between UHECR arrival directions and possible sources 
such as (specific subsets of) active galactic nuclei. In particular a selected 
sample of 39 AGASA events (with E > 4.8 x 10^° GcV) and 26 Yakutsk 
events (with E > 2.4 x 10^° GeV) are claimed to have significant alignments 
within 2.5° with 22 selected BL Lacartae objects (having redshift z > 0.1 or 
unknown, magnitude m < 18 and 6 cm radio flux Fq > 01.7 Jy). The value 
of -P^fn^ is 4 X 10~^ and the penalty factor for making cuts on the BL Lac 
catalogue is estimated to be only 15, yielding a chance probability of 6x 10~^ 
[12]. However the Yakutsk experiment has an angular resolution worse than 
4° at the low (sub-GZK) energy cut made to maximise the coincidences, so 
these cannot be physically meaningful; dropping these events increases the 
chance probability for the remaining coincidences (with AGASA events) 
to 0.15% [13]. Also the assumption that the "energies of the events are 
not important for correlations at small angles" [12] is clearly wrong since a 
mild drop in the energy threshold for AGASA events to £^ > 4 x 10^° GeV 
decreases the significance further by a factor of 5 [13]. Relaxing the cuts 
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Fig. 3. Autocorrelation scan of 57 AGASA events with E > 4 x 10^" GeV shown 
in 4 views [10]. The bottom panel shows the chance probabilities for the observed 
clustering separately for the 30 and 27 events, observed before and after Oct 1995. 



on the BL Lac catalogue yields only 2 coincidences between BL Lacs and 
AGASA 'doublets', with a chance probability of 6.3% [13]. Moreover an 
independent sample of 27 Haverah Park and 6 Volcano Ranch events with 
E > 4 X 10^° GeV do not coincide at all with the chosen 22 BL Lacs [14]. 
Thus, as shown in Fig. 4, there is no justification for the claim that "BL 
Lacertae are sources of the observed ultra-high energy cosmic rays" [12]. 

Finally, the inferred composition of UHECRs (as well as their energies) 
is sensitive to the UHE interaction model used [15] and this realization 
has muddled the picture suggested earlier by Fly's Eye and HiRes of a 
change from heavy nucleus domination at 3x 10^ GeV to nucleon domination 
at 10^'' GeV. Present data is consistent with a mixed composition at all 



Fig. 4. The sky distribution of 57 UHECRs (circles) with E > 4x 10^" GeV observed 
by AGASA, with the 5 'doublets' and 1 'triplet' highlighted. The left panel shows 
the 22 BL Lacs (dots) satisfying specific cuts on redshift, magnitude and 6 cm radio 
flux [12], while the right panel shows all 306 BL Lacs in the catalogue [13]. 

energies. Reanalysis of 'horizontal showers' at Haverah Park requires that 
no more than 50% of UHECRs above 4 x 10^'' GeV can be photons [16] and a 
similar (but weaker) limit has been set by AGASA on the basis that photon- 
initiated showers tend to be muon-poor [17]. However the showering of UHE 
photons is rather complex (because of pair conversion in the geomagnetic 
field and the LPM effect), e.g. it cannot be excluded that the highest energy 
Fly's Eye event (£■ ~ 3 x 10^^ GeV) was in fact initiated by a photon [18]. 

3. Conventional explanations 

No astrophysical object has ever been definitively identified (e.g. through 
7-ray or z/ emission) as an accelerator of high energy nucleons. However 
there do exist energetic objects such as 7-ray bursts, active galactic nu- 
clei, or extended lobes of radio galaxies which satisfy the 'Hillas criterion' 
of being big enough and/or having sufficiently large (estimated) magnetic 
fields to be able to confine UHECRs upto ~ GeV [1]. The issue of 
whether particles can actually be accelerated to the observed energies in 
such objects is an open one, in particular considerations of radiative energy 
losses during the acceleration process set very strong constraints [19, 20]. 
Ultrarelativistic bulk flows in 7-ray bursts (GRBs) remain a possibility, and 
this coupled with the HiRes claim of a GZK cutoff in the energy spectrum, 
have revitalized the suggestion that such cosmologically distant objects are 
the sources of UHECRs [21, 22]. However the energetic requirements are 
formidable — such sources must inject ~ 10^^ erg/Mpc^/yr in UHECRs of 
energy lO^'' — 10^^ GeV in order to match the observed flux, while GRBs 
are observed to emit roughly this much power at far smaller energies of 
0(1) MeV. It is hard to conceive of a particle acceleration mechanism that 
can generate comparable amounts of power in such widely separated energy 
regions. Nevertheless if there is indeed a GZK cutoff, this is probably the 
most attractive possibility for the sources of UHECRs. 
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4. Models involving new physics 

However if the UHECR spectrum does extend without a GZK cutoff 
as indicated by AGASA, then new physics would appear to be required. 
We discuss below only those possibilities which have not been ruled out 
already. (In particular wc do not consider whether the primaries might be 
neutrinos having enhanced interaction cross-sections through new physics 
such as TeV-scale extra dimensions — this cannot explain the observed 
UHECRs, although studies of airshowers can in principle probe such new 
physics [2]. We also disregard the hypothesis that the primaries are new 
light stable hadrons (for which the GZK cutoff energy can be higher) , since 
such particles are excluded by laboratory experiments [23].) 

4-1- Violation of Lorentz invariance 

It was noted long ago that a small modification of the relation between 
momentum and energy in special relativity may undo the GZK cutoff [24]. 
If the clustering of events claimed by AGASA is substantiated by the forth- 
coming high statistics data from Auger, then point sources of UHECRs 
would be implicated. In that case it may be appropriate to invoke such 
violation of Lorentz invariance to explain the absence of the GZK cutoff in 
the spectrum. However it is hard to see how this possibility can be falsified 
since the required violation is so small that it need not manifest itself in any 
other astrophysical or laboratory phenomenon [25]. 

4.2. Z-hursts 

Since at least one species of relic neutrinos must have a mass ^0.1 eV, it 
is attractive to suppose that these provide a target for UHE neutrinos from 
distant sources to annihilate on. This would create 'Z-bursts' with an energy 
of m\/2mi, ~ 4 X 10^'^{mv/l eV)^^ GeV, i.e. in the right energy range to 
be a source for UHECRs [26, 27]. The energy spectrum of the nucleons and 
7-rays resulting from Z decays is well known so a detailed comparison can 
be made with the AGASA data. Agreement is obtained for a relic neutrino 
mass of 0.08 — 1.3 eV [28], however the required UHE u flux is very large, 
taking into account that such light relic neutrinos cannot cluster significantly 
[29]. Although there are no direct bounds yet at these energies, it seems 
implausible that the required extragalactic sources of ~ 10^^ GeV neutrinos 
can exist, given the restrictive bounds on deeply penetrating air showers 
from experiments such as Fly's Eye, as well as AGASA and RICE [30]. 
Moreover it is not clear how such high energy neutrinos can be created in the 
hypothetical cosmic sources — usually this would require the acceleration 
of even higher energy nucleons! 



8 



4-3. Decaying supermassive dark matter 

The possible existence of relic metastable massive particles whose decays 
can create high energy cosmic rays and neutrinos had been discussed [31, 32] 
before the detection of the famous Fly's Eye event with E' ~ 3 x 10^*^ GeV 
which focussed attention on the possible absence of the GZK cutoff. Sub- 
sequently this idea was revived and it was noted that such particles, be- 
ing cold dark matter, would naturally have a overdensity by a factor of 
~ 10^ in the halo of our Galaxy [33, 34]. Hence if their slow decays gener- 
ate UHECRs, all propagation effects will be unimportant (except possibly 
for photons) in determining the observed spectrum and cosmposition, and, 
moreover, there should be a detectable anisotropy in the arrival directions, 
given our asymmetric position in the Galaxy. In order to account for the 
highest energy events with the observed rates, the particle mass should ex- 
ceed mx ^ 10^^ GeV while to match the UHECR flux its lifetime must 
be rx — 3 X 10^^ yr, where (^x is the fraction of the halo dark matter 
in the form of such particles. It is conceivable that such particles exist 
in the 'hidden sector' of string/M-theory [35] and can be produced with a 
cosmologically interesting abundance at the end of inflation [36]. 

The spectra of the decay products is essentially determined by the 
physics of QCD fragmentation so can be calculated e.g. by evolving the 
fragmentation functions measured in decay upto the energies of interest 
using the DGLAP equations [37, 38, 39]. Recent work [40] has included 
electroweak corrections and a more careful treatment of the possible effects 
of supersymmetry. As seen in Fig. 5, the evolved spectrum of nucleons is 
in good agreement with the 'flat' component of UHECR extending beyond 
Eqzk- However the decay photons, which have a similar spectral shape 
are more abundant by a factor of ^ 2, so this model (as well as a similar 
subsequent proposal involving annihilations rather than decays [41]) would 
seem to be ruled out by the experimental bounds [16, 17] on the photon 
component of UHECRs. It is possible that the attenuation of such UHE 
photons in the halo of the Galaxy (through pair production on ^^MHz fre- 
quency radio photons) has been underestimated — the radio background 
intensity assumed for such calculations is very uncertain [42]. However, such 
attenuation would in turn generate a background of low energy 7-rays and 
it is necessary to check that this docs not exceed observational limits, par- 
ticularly from EGRET at ~ 100 MeV. As shown in Fig. 6, such constraints 
can indeed be satisfied although the required increase in the intensity of the 
radio background is rather large [43] . Perhaps one should wait for Auger to 
establish whether photons are indeed ruled out as the UHECRs. 

Detailed calculations have also been made of the expected anisotropy, 
adopting different possible models of the dark matter halo (cusped, isother- 
mal, triaxial and tilted) [44]. The amplitude of the anisotropy is controlled 
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Fig. 5. Fit to the UHECR spectrum beyond Eqzk (excluding the HiRes data) with 
a decaying dark matter particle mass of 5 x 10^^ GeV (dashed line); the dotted line 
is the extension of the spectrum observed at lower energies [39] . 
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Fig. 6. The expected processing of the UHE photon spectrum from decaying dark 
matter in the Galactic halo is shown (left panel) for a radio background 10 times 
higher than is usually assumed and the 7/p ratio for various assumed values of the 
radio background is compared with experimental limits (right panel) [43] . 



by the extent of the halo, while the phase is controlled by its shape. As 
seen in Fig. 7, the amplitude of the first harmonic is ~ 0.5 for a cusped 
'NFW halo but falls to ~ 0.3 for an isothermal halo which is more favoured 
by observations [45], while the maximum is in the direction of the Galactic 
Centre with deviations of up to 30*^ for triaxial and tilted haloes. To detect 
the predicted anisotropy (for the likely case of an isothermal halo with large 
core radius) will require detection of ~ 500 UHECRs by Auger; this should 
also suffice to determine the angular phase to within ±20*^ [44]. 
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Fig. 7. Contour plots of the UHECR sky for (a) cusped, (b) isothermal, (c) triaxial 
and (d) tilted models of the dark matter halo of our Galaxy [44] . 

Another signature of this hypothesis is the expected UHE neutrino flux, 
which was indeed the first to be studied [32]. Recent detailed calculations 
[46] indicate that the expected flux should yield ~ 10 — 40 events/yr with 
E > 10^ GeV in IceCube [49] and a similar number of events in RICE [50]. 
Of course other proposed models for UHECRs also predict associated fluxes 
of UHE neutrinos. However an unique test is the expected flux of neutralinos 
if these are the lightest stable supersymmetric particles [47]. In principle 
EUSO [51] can discriminate such events from neutrino-induced ones [48]. 

5. Conclusions 

After a long hiatus, high energy cosmic rays have again become very 
interesting for particle physicists looking for evidence of physics beyond the 
Standard Model. The source of the highest energy particles in Nature is an 
equally interesting enigma for astrophysicists. As Lemaitre first suggested, 
the origin of such particles may even be linked to the early universe, al- 
though not quite as he imagined! Presently the data are tantalising but not 
sufficient in either quantity or quality to distinguish definitively between 
proposed models. The good news that this will soon be remedied by the 
Pierre Auger Observatory [6]. About 100 water Cerenkov detectors and 
2 fluoresence detectors are now operational, covering an area twice that 
of AGASA, and the full array should be operational by the end of 2005. 
Moreover the ambitious space-based experiment EUSO [51], scheduled for a 
3-year engineering flight on the International Space Station 'Columbus' in 
2007, will provide another substantial increase in collecting power. This is 
an exciting time for cosmic ray physics and we look forward to the surprises 
that Nature has in store for us. 
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